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Abstract 



We study the pion form factor F 71 " 77 * (Q 2 ) in the light -cone sum rule approach, 
accounting for radiative corrections and higher twist effects. Comparing the results 
to the CLEO experimental data on i* 171 " 77 (Q 2 ), we extract the the pion distribution 
amplitude of twist-2. The deviation of the distribution amplitude from the asymp- 
totic one is small and is estimated to be 0*2(11) = 0.12 ± 0.03 at u = 2.4 GeV, in the 
model with one non-asymptotic term. 

The ansatz with two non-asymptotic terms gives some region of ai and 04 , which 
is consistent with the asymptotic distribution amplitude, but does not agree with 
some old models. 
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1 Introduction 



The production process of one neutral pion by two virtual photons, 7*7* — > tt°, plays 
a crucial role in the studies of exclusive processes in quantum chromodynamics. Being 
one of the simplest exclusive processes, it involves only one hadron and relates directly 
to the pion distribution amplitude ffl]. At large photon virtualities, we can calculate the 
form factor using perturbative QCD and obtain important information on the shape of 
the pion distribution amplitude from the experimental data. 

In general, the pion distribution amplitudes serve as input in the QCD sum rule 
method and allow the calculation of many form factors (for example heavy-to-light form 
factors B — > 71 and D — > 71) and hadronic coupling constants (for example gs*B-K and 
9d*Dt^)- We refer the reader to reviews || |7| and recent studies ||. 

Recently, the CLEO collaboration has measured the 77* — ► tt° form factor. In this 
experiment, one of the photons is nearly on-shell and the other one is highly off-shell, with 
a virtuality in the range 1.5 GeV 2 - 9.2 GeV 2 f§ ft We study the possibility of extracting 
the twist-2 pion distribution amplitude from the CLEO data. 

The pion-photon transition has been the subject of many studies in framework per- 
turbative QCD |13[], lattice calculations |18j and QCD sum rule method |12|, [L4|, [13], [19 . 



In the present paper, we use the light-cone sum rule (LCSR) method for calculating 
the form factor of the process 77* — > tt° . The method of LCSR has been suggested in 
H |l(J and consists of the operator product expansion (OPE) on the light-cone @, [lj] 



combined with the QCD sum rule technique |L2]. 

The first attempt to calculate this form factor by using the LCSR method has been 



reported in ||19|| . The twist-2 and twist-4 contributions were calculated to leading order 
without accounting for perturbative QCD effects. The radiative QCD effects are usually 
large (20%) and allow to fix the normalization scale dependence of involved parameters. 

In this paper we analyze the LCSR for the form factor 77* — > 7r° at next-to-leading 
order in as (NLO). We derive 0(ots) radiative corrections to the spectral density of 
the twist-2 operator. We combine the twist-2 contribution at NLO with higher twist 
contributions (twist-4) in order to analyze the LCSR for the form factor of the process 
77* — > it numerically. 

Using CLEO experimental data, we extract the parameters of the twist-2 distribution 
amplitude. The distribution amplitude appears to be very close to the asymptotic distri- 
bution amplitude. The deviation of the distribution amplitude from the asymptotic one 
is small and is estimated to be 



a 2 (n) = 0.12 ±0.03 at // = 2.4 GeV, 

using pion distribution amplitude with one non-asymptotic term. This result on a 2 agrees 
well with recent analysis on electromagnetic pion form factor p6| . 

The pion distribution amplitude with two non-asymptotic terms is also extracted, 
giving some region in the a 2 , c^-plane. Theoretical systematical uncertainties are dominant 
over experimental statistical ones and systematical uncertainties are strongly correlated, 
defining the allowed parameter space of 02,04. The extracted region for the distribution 

1 There exist also older results by the CELLO collaboration [El. 
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amplitude is in qualitative agreement with results derived in [15L ITO, [L7|, E6| , where it has 



been claimed that the pion distribution amplitude is very close to the asymptotic form. 
However, our region does not overlap with the pion distribution amplitudes suggested in 
|nj and §. ' 

The paper is organized as follows. In the next section we discuss the general framework 
of calculating the 7177* form factor with LCSR. In section 3 we calculate a spectral 
density at LO and NLO and present the final sum rule at NLO. In section 4 we perform a 
numerical analysis and discuss a procedure to extract the parameters a 2 and 04 including 
the estimation of the systematic and statistical uncertainties. 



2 The method of calculation 

We start with the correlator of two vector currents j^/ v = (^uj^/ u u — \d^^/ u d): 

J d A xe~^ x < 7r°(0)|T{j M (x)j i ,(0)}|0 >= ie^ af} q?q$F^*( qi ,q 2 ), (1) 

where (71,(72 are the momenta of the photons. If both virtualities Si = q\ and s 2 = q\ 
are large and Euclidean, — si,— s% Aq C d, the correlator can be expanded near the 
light-cone (x 2 — > 0). The leading twist-2 contribution to the correlator (|l|) is [p]] 

3 J s 2 (u — 1) — SiK 

where (p n (u) is the pion distribution amplitude of twist-2 defined through 

1 

7^75 4(0)10) = -iq^U fdu V Me luqx . (3) 



In principle, the higher twist contributions can be calculated using the OPE on the light 
cone. 

However, in the CLEO experimental data, one of the virtualities is small, i.e. s 2 — > 0. 
A straightforward OPE calculation is not possible and we have to use analyticity and 
duality arguments. 

Since the form factor F wl 7 (s\, s 2 ) is an analytical function in both variables, we can 
write the form factor as a dispersion relation in s 2 : 



mi - s 2 J s - s 2 

f S 



(4) 



The physical ground states are p and u> vector mesons. We use the zero-width approx- 
imation and define the matrix elements of electromagnetic currents, assuming isospin 
symmetry: m p ~ m u ; |(7r°(p)|j M |^(g 2 )) - (^°(p)\j fl \p°(q2)) = ^ e^af} e u q% q% F p7r (si); 

3 (u\j u \0) ~ (p°|ji/|0) = ^m p e* v . Here e v is the polarization vector of the p meson; f p is 
its decay constant. 
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The spectral density of higher energy states p h (si, s) is derived from the QCD-calculated 
expression for Fqqq (si, s) by using usual semi-local quark-hadron duality for s > s 

Jds(any function) / s) = J ds(any function)! Imi^SV, s). (5) 

We equate the dispersion relation (|j) with the QCD expression at large s 2 . Using the 
dispersion relation for the QCD function Fqq^ (si, s 2 ), we obtain 

V2f p F^( Sl ) = 1 ImFgljs^s) 
m 2 p - s 2 TT J s - s 2 

The next step is to perform a Borel transformation in s 2 . We finally get the LCSR for 
the form factor F p7r (si): 

V2f p F^( Sl ) = - JdslmF^(s 1 ,s)e : ^ 1 , (7) 
o 

where M is a Borel parameter. 

Substituting (|7]) and the duality approximation (|5|) into and taking the s 2 — > 
limit we obtain the sum rule for the form factor F 77 * 7r (si) | I9| : 

F^ ( Sl ) = — da1mF%2 (s,, s) + - /-ImF^ ( Sl ,s). (8) 

p so 

We will use this expression as our basic sum rule for the numerical analysis. 



3 Born term and QCD radiative correction 

The next step is to calculate the spectral density at LO and NLO. 

Calculating twist-2 contributions to F 71 " 7 * 7 * (si, s 2 ), only one distribution amplitude en- 
ters, the pion distribution amplitude, ip^(u), defined by ([$]). As a result, the F 7r7 * 7 *(si, s 2 ) 
can be written as the convolution of the hard amplitude T(s\, s 2 , u) and the distribution 
amplitude (pn(u): 

i 

pn-f-f ( Sl; S2 ) = Jdu<pn(u) T(si, s 2 ,u). (9) 



The hard amplitude - it plays the role of the Wilson coefficient in OPE - is calcula- 
ble within perturbative theory, the pion distribution amplitude <p n (u) contains the long- 
distance effects. 

The theoretical spectral density at s 2 > and si < can be calculated from 

1 } 1 

Pqcd(si,s 2 ) = -lm S2 F nl 7 (s 1 ,s 2 ) = f w du<p„(u) -Im S2 T(si, s 2 , u). (10) 

TT J TT 
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The contribution of twist- 2 at LO approximation (0) M is 



^^1^) = ^ 1 iu and P L Q °cd = ^Z^- ) Muo). (11) 



where uq 



S2-S1 



The 0(a s ) correlation function for 7*7* — > 7r° was considered in |20| (see also [[H], ^2 



Here we use this result and present it in a form which is useful for our further calculations 



c 



a s (M)C F 
3(s 2 -si) 2 



7T / ' 



J QOD - J QCD + 1 QCD > U^J 



Toot ( s i> S2,w) = c(a L + ai Li + 02-^2) , with L n = lQg — — — . (13) 

^ U — Uq 

The coefficients ao are expanded in terms of log n (— no): 

«o = b 1 + 61 /1 + 6 2 h with i n = log n (-u ). (14) 

The coefficients are: 

-S^M + S 2 Si(l +U) 

°2 - ^—7 , (.loj 

2-u(s 2 - 

h = -2& 2 (L M + 1)- 52 , 

2(s 2 -si)m 

60 = L^ + Lj- \ 2 (L, + l)b] -§(3 + L M ), 



a 2 



■§2 Si 



2 2(s 2 -si) 2 u(m-1)' 



s2(l - u) 2 + s 2 u 2 + s 2 si (-3 + 2u - 2u 2 
ai = — 2a 2 L M H 



with 



2(s 2 - Si ) 2 (m- l)u 



log( ); L u = log(l - m ); 6- 



s 2 - si " ' ' 2 (s 2 - si) 2 (w - 1) 

Now we take the imaginary part of this expression in the energy region s\ < and 
S2 > 0. The nontrivial imaginary part comes from the functions L n and l n . We collect all 
useful formulae in Appendix A. The combined result is 

- ImT™( Sli s u) = -J F (—LJ\ ( 6l + 2 log(n ) 6a) (16) 

7T V y \U — UqJ 

- 8{u- u ) (b + log(w )(6i + at) + log 2 (w )(&2 + a 2 ) - 7r 2 + ^2)) 
n , x / / 1 \ . f2log\u -u\\ \ \ 



U — Uq / + \ U — Uq 
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+/ 



where we use the () + operation, which is defined inside the integral as: 

J duf{u)(g(u,u f) + = J du (f(u) - f(uo))g{u,u ). 

We have checked that the dispersion integral with the imaginary part (|I6|) reproduces the 
hard amplitude flTB"!). 



The expression (1^) is universal and could be used with any distribution amplitude. 
In the present paper we expand the distribution amplitude in Gegenbauer polynomials, 
keeping the first three terms (see also next section for details): 

ri w2 (u,fi) = 6w(l - u) (l + a 2 (n) Cl /2 (2u- 1) + a 4 (/x) Cl' 2 (2u - 1) + . . .) . (17) 
The direct integration over variable u gives the theoretical spectral density 

P$g(« a , = ^ a ' ( 2 f F {M*2, *i) (18) 

+ a 2 (/i) A 2 (s 2 , si, fi) + a 4 (/i) A 4 (s 2 , s 1: /x) 
with the coefficients 

A = . S2S \ 3 (-15 + 7 r 2 -31og(-^)\ (19) 

(s 2 - Sl y \ s x j 

A 2 = 82 , f-25s 2 3 -8 (95 + 3tt 2 ) s 2 2 Sl -36 (25 + 2tt 2 ) s 2Sl 2 

4(s 2 -s a ) V V J V J 

- 12 (5 + 2 7T 2 ) si 3 + 12 si (s 2 2 + 3 s 2 Sl + si 2 ) ^—25 log(^-) + 6 log(-^) 2 

A A = 7 f-91s 2 5 - 2 (541 3 + 75tt 2 ) s 2 4 si- 125 (541 + 12tt 2 ) s 2 3 s^ 

10(s 2 -si) V v ' v ' 



- 100 (901 + 30 7T Z J s 2 z - 150 (193 + 10 tt z J s 2 - 15 ( 109 + 10 tt z ) Sl " 
+ 30 si f s 2 4 + 10 s 2 3 si + 20 s 2 2 si 2 + 10 s 2 si 3 + Si 4 ) x 



^_911 g ( ^!) + 151 0g( _|) 2 ^. 



As a result the asymptotic contribution of twist-2 operators at NLO has a very simple 
form 

PQCD = ~ rg-(H : rz (-15+7T -31og( ) JJ (20) 

(S 2 - Si) V 127T v Si // 

We note here that the spectral density contains double logarithms a s log(— ^-) 2 . Numer- 
ically, these are moderate in our LCSR with s 2 « M 2 « s = 0(1 GeV 2 ) and -s x < 10 
GeV 2 (as in the case of CLEO data). For higher virtualities, —Si ^> 10 GeV 2 , the resum- 
mation of the double logarithms will be necessary. 

Finally, we obtain LCSR by combining the general expression (|j) with the formulae 
for the twist-2 spectral density at NLO derived in this paper (|H], [IB], [B5]) with the twist-4 



contribution taken from the literature 19 
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4 Numerical results 



We use the following parameters from the Particle Data Group [23| for the numerical 
analysis: f w = 132 MeV, f p = 216 MeV and m p = 770 MeV |23[. For the running of 
the QCD coupling constant as, we use the two-loop expression with N f = 3 and A^ 4 ) = 



380 MeV which corresponds to as(Mz) = 0.118 |23|, after matching twice the QCD 
coupling constant at the quark-antiquark thresholds /i c5 = 2.4 GeV and n b i = 10 GeV. 

The distribution amplitude can be expanded in terms of Gegenbauer polynomials 
^n( u ) — 6w(l — u)C^ 2 (2u — 1). Arguments based on conformal spin expansion [|Il[] allow 
us to neglect higher terms in this expansion. We adopt here the ansatz which consists of 
three terms, assuming that the terms a„>4 are small 

<p v (u, Ho) = * («) + a 2 (fx )^ 2 (u) + a 4 ()"o)*4(«)- (21) 

The asymptotic distribution amplitude ip n (u) = 6w(l— w) is unambiguously fixed ||. The 
terms n > describe non-asymptotic corrections. 

We may consider two different approaches to confront the sum rule with the experi- 
mental data. The first possibility is to use the existing values for 02, 04 JIT], H in our sum 
rule and to compare the calculated form factor with the CLEO experimental data. The 
second approach is to treat 02, 04 as unknown parameters and to extract them from the 
CLEO data. We shall discuss both possibilities, starting with the first one. 



4.1 The asymptotic, BF- and CZ- distribution amplitudes 

Previous extractions of the parameters a 2 , 04 can be divided into three classes: 

1) Chernyak and Zhitnitsky (CZ) obtained the coefficients a 2 (fi ) =2/3 and 04(^0) = 
[||] at the scale /x = 0.5 GeV. 

2) Braun and Filyanov (BF) extracted the coefficients a 2 (/io) — 0.44 and a 4 (/i ) = 0.25 
at the scale /i = 1 GeV. 



3) Some groups argued that the wave function is very close to the asymptotic one, the 



coefficients are very close to a 2 = and a 4 = (see for example 0, [13|, |H| [16], fT?] ) 
The QCD evaluation of a 2 and 04 at different scales /i gives 





0.7 GeV 


1.0 GeV 


1.5 GeV 


2.4 GeV 




a 2 


04 


a 2 


04 


a 2 


(24 


a 2 


Cl4 


BF 


1 1 1 






0.44 


0.25 


0.33 


0.17 


0.28 


0.13 


CZ 




2/3 









0.33 


-0.008 


0.28 


-0.009 



Now we calculate the form factor using these values. We use the normalization point 
fi = 2.4 GeV, which corresponds to the virtuality of the photons in the central region of 
the experimental data. This scale is also used in many B — > it calculations |§ . The lower 
scale h = 1.5 GeV is taken in order to check the sensitivity of the results to the variation 

of /!. 
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The sum rule depends on the Borel mass M and the threshold energy sq. The depen- 
dence on both parameters is small. The variation of so by 20% gives deviations in the 
form factor of less than 2%. Fig. [I] shows the Borel mass dependence of the form factor at 
various Q 2 (Q 2 = — si). We note that the dependence is very small for virtualities around 
Q 2 2 -f- 3GeV 2 , where the experimental data is concentrated. For other Q 2 the Borel 
dependence is still moderate (±2 4%), showing a good quality of the sum rule. We use 
in the calculations s = 1.5 GeV 2 and M 2 = 0.7 ± 0.2 GeV. 

Now we are ready to compute the form factor using BF-, CZ- and asymptotic dis- 
tribution amplitudes. In Fig. |2| we show the contributions of the asymptotic and non- 
asymptotic parts of the distribution amplitude at // = 1.5 and 2.4 GeV; the parameters 
a 2 and CI4 are normalized to 1 in the whole region of Q 2 . The twist-4 parameter 5 2 (fi) is 
fixed at fi = 1 GeV to be 0.2 and scaled by the renormalization group equation with the 
one-loop anomalous dimension. We see that the asymptotic contribution is too small in 
order to fit experimental data. 

The results with BF- and CZ- distribution amplitudes are presented in Fig. [| and are 
too large to describe CLEO data. 



4.2 Extraction of the distribution amplitude from CLEO data 

We use a numerical nonlinear fit procedure to estimate best values for Gegenbauer coef- 
ficients. 

First, we study the ansatz for the pion distribution amplitude with only one non- 
asymptotic term a 2 . Using this ansatz and comparing our results with CLEO results we 
obtain 

a 2 = 0.12 ±0.03 at /i = 2.4 GeV. (22) 

This result shows that the pion distribution amplitude is very close to the asymptotic 
form. This conclusion is in agreement with recent analysis on electromagnetic pion form 



factor [pg] an d also with results based on QCD sum rules presented in |L5|, [16], |17 |. 

Let us comment on the normalization scale dependence of our result, which enters 
through the QCD coupling constant and logarithms in the radiative correction. For the 
extraction of a 2 (/i), we may choose any reasonable normalization point in the interval 
1 < /i < 3 GeV. Different input values of ji will give different pairs of a 2 (/i). We have 
checked that all extracted values of a 2 are in agreement with the renormalization group 
equation. The details about the running of a n (/x) with fj, are discussed in Appendix B. 

We have also studied the ansatz for the pion distribution amplitude with two non- 
asymptotic terms a 2 and 04. Using this ansatz and comparing our results with CLEO 
results we obtain: a 2 = 0.19 and 04 = —0.14 at fi = 2.4 GeV as the fit parameters. 
The form factor calculated with the central values of the extracted parameters 02(2.4) 
and a 4 (2.4) is shown in Fig. 

Now let us discuss systematic uncertainties. There are many sources of systematic 
uncertainties. One is the high order QCD perturbative corrections. Taking into account 
the size of one loop QCD correction (20%), we estimate two loop (and higher order 
correction) to be of order 0.2 • 0.2 = 0.04(4%). Additionally, there are power corrections 
of twist higher than 4. We assume them to be 20% of the twist-4 contribution, which itself 
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contributes about 20 — 30% to the twist-2 term. The effect of the twist-4 contribution is 
shown in Fig. [5]. We also adopt zero width approximation for p and u. The uncertainty 
induced by this assumption is small, a few percent^, since the form factor F p7T does 
not appear in the final sum rule. The duality assumption introduces another source of 
uncertainties which is difficult to estimate. We assume that the variation of M 2 and So 
gives us a rough idea about the size. 

To fix statistical uncertainties we apply the fit procedure on a statistical set of data 
tables, where the experimental points are randomly displaced within the given errors 0. 

Our results are presented in Fig. |G], where we show a parameter space of the twist-2 
distribution amplitude, which is obtained by comparing CLEO results with our sum rule. 
Besides of statistical uncertainties, we account for the possible theoretical uncertainties. 
We combine theoretical uncertainties together and write the theoretical form factor as 
(with 95% CL): 



"T7T77 



= Flr 2 (si)(l ± 5%) + FXrXs^l ± 20%) + at F a2 ( Sl ) + a{ u F a4 ( Sl ). (23) 



The F^~ 2 ^ denote the asymptotic contribution of twist-2/twist-4 distribution ampli- 
tudes, F a2 ,F a4 are the normalized contributions of the higher Gegenbauer polynomials. 

The width of the allowed region ( Fig. ||) is due to experimental-statistical uncertain- 
ties, whereas the length is due to theoretical-systematical ones. 

As we see form Fig. ||, theoretical-systematical uncertainties are dominant over exper- 
imental statistical ones, giving a correlation between the parameters 04. 

Fig. [7] shows a more conservative scenario where we assume an uncertainty of ±8% in 
the twist-2 term and again ±20% for the twist-4 contribution. 

It is worth to mentioned that the allowed space of a 2 and 04 does not overlap with 
Braun-Filyanov and Chernyak-Zhitnitsky distribution amplitudes [|[ [LI]. These distri- 



bution amplitudes are beyond the 95% CL level region of the allowed parameter space 
for a,2, 04. At the same time our region for CL2, CI4. IS consistent with results presented in 
15|, ITBl [IT], |26|, where it has been claimed that the pion distribution amplitude is very 



close to the asymptotic form. 

The region presented in Fig. [7j does not favor any central value. But in order to 
quantify the region, we may roughly identify following values for a 2 and a 4 : a 2 = 0.19 ± 
0.04(stat.) ± 0.09(syst.) , a 4 = -0.14 ± 0.03(stat.) =f 0.09(syst.) 

We also observe that F 71 " 77 * (Q 2 ) form factor is more sensitive to the sum of 02 and 04, 
than to the 

O2 — ^4- 

Really, we found that only one linear combination of 0,2 and 04 is 
determined well, it is a 2 + 0.604 = 0.11 ± 0.03 at /i = 2.4 GeV. 

It would be very useful to study some observable, which is sensitive to 

O2 — O'Ai 

for 

example, the coupling constant g pu}1T . LCSR for this coupling is related to ^(1/2), which 
is sensitive to the combination d 2 — d\- 

In general, the result for the model with one non-asymptotic term (^) is very reliable. 
At the same time, we have very useful constraint on the model with two non-asymptotic 
terms. Although, an additional constraint on a 2 — CI4 would be extremely useful. 

2 In fact, at the — s\ > 4 GeV 2 the resonance part in (8) (dual to the form factor F p7T ) contributes less 
then 20 — 30%. The finite width of p meson gives a deviation of order 10% ]19| , therefore, we are left 
with 2 — 3% uncertainty in the region — s\ > 4 GeV 2 . It is worth to note that this energy region gives 
the same values for 02, 04 as the overall sample of data. 
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5 Conclusions 



We have studied the pion form factor F 71 " 77 * (Q 2 ) in the light-cone sum rule approach, 
including radiative corrections and higher twist effects. Comparing the results to the 
CLEO experimental data on F 71 " 77 * (Q 2 ), we have extracted the pion distribution amplitude 
of twist-2. The deviation of the distribution amplitude from the asymptotic one is small 
and is estimated to be a 2 (/i) = 0.12 ± 0.03 at // = 2.4 GeV, in the model with one 
non-asymptotic term. This result on a 2 is in agreement with very recent analysis on 
electromagnetic pion form factor | 26| . The ansatz with two non-asymptotic terms gives 



some region of a 2 and 04, which is in qualitative agreement with results derived in 
15], O, 53], but does not agree with CZ and BF models [H, 111 . 



Since the theoretical systematical uncertainties are dominant over experimental sta- 
tistical ones it will be very useful to estimate the power corrections and the higher order 
radiative corrections to the correlator in future. Also, it would be very useful to calculate 
LCSR for some observable, which is sensitive to 02 — 04, for example, the coupling constant 
gpun, and to confront LCSR results with experiment. 
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6 Appendix A 



To calculate the imaginary part, we use following identities: 
— IhiLq = 5(u — Mo); 

7T 



ReL, 



1 



o 



U — Uq 



-lmL 1 = -<d(u -u)( — - — ) + S(u - u ) log \u \; 

7T \U — Uq/ + 

-lmL 2 = -Q(u - u) ( 21og ^° — —\ - 5(u - w ) ( ?r - log 2 |w | ] ; 

7T \ « - «0 / + V 3 / 

-Imh = -9(«o); 

7T 

Re/i = log|«o|; 
-ImZ 2 = -2©(« ) log(wo); 

7T 

Re/ 2 = log 2 \uq\ — 7t 2 9(m ). 

7 Appendix B 

The distribution amplitude <p n can be expanded in terms of Gegenbauer polynomials 
^n( u ) — 6m(1 — u)C% 2 (2u — 1). In NLO, the evolution of the distribution amplitude is 
given by |22|: 

/i) = E «»(tt>) exp f - / da^A U n (u) + ^ £ d k MMu)) (24) 

with a = 1. The coefficients c?^(/i) are due to mixing effects, induced by the fact that the 
polynomials ^ n {u) are the eigenfunctions of the LO, but not of the NLO evolution kernel. 
The QCD beta-function f3 p3| and the anomalous dimension 7™ of the n-th moment a n (n) 
of the distribution amplitude have to be taken in NLO. Explicitly at NLO, the exponent 
in (M) is 



,,n _ n 



Ait 



The anomalous dimensions E5 



are 

•^-s^Cs)'* (26) 

with 

7o° = 0, 7i° = 0, 

830 

-JV F . (27) 



7o° 


= o, 


7? 


= , 




100 




34450 


7o 2 


9 ' 


7 2 


243 




728 




662846 


7o 4 




it 




45 ' 


3375 



— N F , 
51 



2025 
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Np is a number of active flavours. The beta-function coefficients are defined in a standard 
|. The NLO mixing coefficients are E2, 04 



way 



rtfc 



7o fe - 7o n - 2/3o 



V 



2/30 



:28i 



/ 



where the numerical values of the first few elements of the matrix M n ^ are 

M 02 = -11.2 + 1.73JV>, M 04 = -1.41 + 0.565iV F , M 24 = -22.0 + 1.65iV F . (29) 
The QCD evaluation of a 2 , a 4 to the scale \i = 2.4 GeV gives 

02(2.4) = 0.28, a 4 (2.4) = 0.13 
for the Braun-Filyanov distribution amplitude and 

02(2.4) = 0.28, 04(2.4) = -0.009 
for the Chernyak-Zhitnitsky distribution amplitude. 



12 



References 

[1] G.P. Lepage and S.J. Brodsky, Phys. Lett. B87 (1979) 359; Phys. Rev. D22 (1980) 
2157. 

[2] A.V. Efremov and A.V. Radyushkin, Phys. Lett. B94 (1980) 245; Teor. Mat. Fiz. 
42 (1980) 147 

[3] V.L. Chernyak and A.R. Zhitnitsky, JETP Lett. 25 (1977) 510; Sov. J. Nucl. Phys. 
31 (1980) 544; Phys. Rep. 112 (1984) 173. 

[4] CLEO collaboration, V. Savinov, Contribution to the Intern. Conference Photon 97, 
Edmond aan Zee, Netherlands, (1997), |hep-ex/ 9707028 ; 



CLEO Collaboration (J. Gronberg et al), preprint CLNS-97-1477 (1997), |he£ 



ex/9707031 



[5] CELLO collaboration, H.-J. Behrend et al Z. Phys. C 49 (1991) 401. 
[6] V. Braun, Preprint NORDITA-98-1P (|hep-ph/9801222j ). 



[7] A. Khodjamirian and R. Riickl, |hep-ph/ 980 1443|, in Heavy Flavours II, eds. A.J. 
Buras and M. Lindner (World Scientific, Singapore, 1998), p. 345. 

[8] A. Khodjamirian, R. Riickl, S. Weinzierl, O. Yakovlev Phys. Lett. B 410 (1997) 275; 
E. Bagan, P. Ball, V. M. Braun, Phys.Lett. B417 (1998) 154-162; 
P. Ball JHEP 9809 (1998) 005, |hep-ph/9802394] ; 

S. Weinzierl, O. Yakovlev, [hep-ph/9712399| , in "Progress in Heavy Quark Physics", 
ed. T. Mannel; 

A. Khodjamirian, R. Riickl, S. Weinzierl, C. Winhart, O. Yakovlev, |hep-ph/00012971 ; 
R. Riickl, A. Schmedding, O. Yakovlev, in preparation; 
A. Schmedding, Diploma Thesis (1999), unpublished. 

[9] I.I. Balitsky, V.M. Braun and A.V. Kolesnichenko, Nucl. Phys. B312 (1989) 509. 

[10] V.L. Chernyak and A.R. Zhitnitsky, Nucl. Phys. B345 (1990) 137. 

[11] V.M. Braun and I.E. Filyanov, Z. Phys. C44 (1989) 157; 
V.M. Braun and I.E. Filyanov, Z. Phys. C48 (1990) 239. 

[12] M.A. Shifman, A.I. Vainshtein and V.L Zakharov, Nucl. Phys. B147 (1979) 385, 448. 

[13] I.V. Musatov, A.V. Radyushkin, Phys. Rev. D56 (1997) 2713. 

[14] A. V. Radyushkin, |hep-ph/9707335| ; 

A. V. Radyushkin and R. Ruskov, |hep-ph/9706518| ; 

A. V. Radyushkin, Few Body Syst. Suppl. 11, 57 (1999) ||hep-ph/9811"225 



A .V. Radyushkin, R.T. Ruskov, Nucl. Phys. B481 (1996) 625. 

[15] S.V. Mikhailov, A.V. Radyushkin, Sov. J. Nucl. Phys. 49 (1989) 494; Sov. J. Nucl. 
Phys. 52 (1990) 697; Phys. Rev. D 45 (1992) 1754. 



13 



[16] A. Bakulev, S. Mikhailov Z. Phys. C 44 (1995) 831; Mod. Phys. Lett. All (1995) 
1611; Phys. Let. B 436 (1998) 351. 



[17; 

[18 
[19 
[20 
[21 
[22 

[23 

[24; 

[25 
[26 



V.M. Belyaev, Mikkel B. Johnson, Mod. Phys. Lett. A13 (1998) 2909. 

D. Daniel et al, Phys. Rev. D 43 (1991) 3715. 
A. Khodjamirian, Eur. Phys. J. C6 (1999) 477. 

E. Braaten, Phys. Rev. D28 (1983) 524. 

F. del Aguila and M.K. Chase, Nucl. Phys. B193 (1981) 517. 

E. P. Kadantseva, S.V. Mikhailov and A.V. Radyushkin, Sov. J. Nucl. Phys. 44 
(1986) 326. 

Particle Data Group, Eur. Phys. J. C 3 (1998) 1-794. 

F. M. Dittes and A.V. Radyushkin, Phys. Lett. B134 (1984) 359; 
M.H. Sarmadi, Phys. Lett. B143 (1984) 471; 

S.V. Mikhailov and A.V. Radyushkin, Nucl. Phys. B254 (1985) 89. 

A. Gonzales- Arroyo, C. Lopez and F.J. Yndurain, Nucl. Phys. B153 (1979) 161. 

V. M. Braun, A. Khodjamirian and M. Maul, Phys. Rev. D61, 073004 (2000), jKep 



ph/9907495 



14 



0.4 



0.5 



0.6 



0.7 



0.8 



0.9 



M 2 (GeV 2 ) 



Figure 1: The form factor Q 2 F 1 ~ / *' !T (Q 2 ) as a function of the Borcl 
parameter M 2 for Q 2 = 1, 2, 9 GeV 2 (from lower to upper line). 
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Figure 2: The contributions of the asymptotic distribution ampli- 
tude (solid line) and non-asymptotic terms of the first two Gegen- 
bauer polynomials with a 2 = 1 in the whole Q 2 -region (dashed 
line) and 04 = 1 (dotted line) to the form factor Q 2 F 77 * 71 (Q 2 ) as 
a function of Q 2 at the normalization point fi = 2.4, 1.5 GeV (up- 
per, lower curves from a bunch of two curves correspondently) The 
experimental data is taken from M . 
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Figure 3: The form factor Q 2 F 11 n (Q 2 ) calculated with differ- 
ent distribution amplitudes: the Braun-Filyanov (dashed lines), 
Chernyak-Zhitnitsky (dotted lines) and our results extracted from 
CLEO data at \x = 2.4 GeV (upper lines from bunches of two lines) 
and fi = 1.5 GeV (lower lines). The experimental data is taken 
from M. 
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Figure 4: The form factor Q 2 F 11 * 1T (Q 2 ) with the distribution 
amplitude extracted from CLEO data at \i = 2.4 GeV (solid lines). 
The contribution of ci2 = 0.19 (dashed line) and 04 = —0.14 (dotted 
line) are the central values of our nonlinear fit. The experimental 
data is taken from Iffl. 
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Figure 5: The form factor Q 2 F 7 ' r * 7r (Q 2 ) as a function of Q 2 at 
best extracted values of a?, and 04 at // = 2.4 GeV with (solid line) 
and without (dashed line) twist-4 contribution. The experimental 
data is taken from ||. 
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Figure 6: The parameter space of (02, 04) pairs extracted from 
3000 randomly chosen sets of data allowed by the experimental 
statistical uncertainties J| as well as by the theoretical systemat- 
ical uncertainties fl2"3"|). Countor-lines show 68% (solid line) and 
95% (dashed line) confidential regions. Bold dots show the param- 
eter pairs for asymptotic (circle) and Chernyak-Zhitnitsky (square) 
distribution amplitudes. 
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Figure 7: The parameter space of (a 2 , a 4 ) pairs extracted from 
3000 randomly chosen sets of data allowed by the experimental 
statistical uncertainties Q as well as by the theoretical systematical 
uncertainties estimated in a very conservative way. Countor-lines 
show 68% (solid line) and 95% (dashed line) confidential regions. 
Bold dots show the parameter pairs for asymptotic (circle) and 
Chernyak-Zhitnitsky (square) distribution amplitudes. 
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